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Abstract 

Cosmic vortons are closed loops of superconducting cosmic strings carrying current 
and chargé. Despite a large number of studies the existence and stability of cosmic 
vorton solutions is still an open problém. Numerical simulations of the nonlinear ŕield 
theory are dimcult to perform in (3+l)-dimensions, due to the existence of multiple 
length and time scales. In this paper we study a (2+l)-dimensional analogue of cosmic 
vortons, which we refer to as kinky vortons, where the cosmic string is replaced by 
a kink string. Many of the expected qualitative aspects of cosmic vortons transfer to 
kinky vortons, with the advantage that several approximations used in the study of 
cosmic vortons can be replaced by exact results. Furthermore, the numerical study of 
kinky vortons requires less computational resources than cosmic vortons, so a number 
of issues can be addressed in some depth. The rádius of the kinky vorton is deter- 
mined as a function of the chargé and winding number, and it is shown that the chiral 
limit is a repulsive fixed point. Stability to both axial and non-axial perturbations is 
demonstrated in the electric and chiral regimes, though surpringly long lived ringing 
modes are observed. Kinky vortons which are too magnetic are shown to suffer from 
a pinching instability, which results in a reduction in the winding number and can 
convert magnetic into electric solutions. 



1 Introduction 



Cosmic strings are topological defects which may have formed during a phase transition 
in the early universe (for a review see [H]). Witten [15] pointed out that if the field of a 
cosmic string is coupled to another complex scalar field then a non-dissipative current can 
form to produce a superconducting string. A cosmic vorton |6j is a stable closed loop of 
superconducting cosmic string which carries both current and chargé. The idea is that the 
current and chargé on the string provide a force to balance the string tension and prevent 
its collapse. 

The cosmological consequences of cosmic vortons have been well studied [2] and it has 
been proposed that they play a role in a number of cosmological phenomena, such as galactic 
magnetic fields, high energy cosmic rays, gamma ray bursts and baryogenesis. However, 
the existence and stability of cosmic vortons as classical field theory solutions is still an 
open problém, and convincing numerical evidence is still lacking. Numerical simulations of 
the relevant nonlinear field theory are difficult to perform in (3+l)-dimensions and results 
are limited, mainly due to the existence of multiple length and time scales. Even in the 
simplest model, which is the global version of Witten's U (ľ) x U(ľ) model [15J, cosmic 
vorton solutions have not yet been constructed. The only field theory computation to date 
[TU] is in a modified version of this global theory, in which the interaction term between 
the two complex scalar fields is replaced by a non-renormalizable interaction, designed to 
make the numerical problém more tractable. Even in this case the vorton solution presented 
eventually decays, though this may be due to numerical issues. 

In this paper we study a (2+l)-dimensional analogue of cosmic vortons, in which the 
cosmic string is replaced by a kink string. For this reason we name these objects kinky 
vortons. We shall show that many of the expected qualitative features of cosmic vortons can 
be demonstrated explicitly for kinky vortons, with the advantage that several approximations 
used in the study of cosmic vortons can be replaced by exact results. Furthermore, the 
numerical study of kinky vortons obviously requires less computational resources than cosmic 
vortons, so a number of issues can be addressed in detail. In particular, the simulations 
presented in this paper reveal some of the computational difficulties involved and help explain 
why similar numerical investigations in (3+l)-dimensions have not yet been successfully 
per for med. 

Our results include the determination of the rádius of a kinky vorton as a function of 
the chargé and winding number, and the demonstration that the chiral limit is a repulsive 
fixed point. Stability to both axial and non-axial perturbations is exhibited in both the 
electric and chiral regimes, though surpringly long lived ringing modes are observed. Kinky 
vortons which are too magnetic are found to suffer from a pinching instability, which results 
in a reduction in the winding number allowing conversion of magnetic into electric solutions. 
Finally, we discuss the implications of our results on kinky vortons to possible future work 
on cosmic vortons. 
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2 The model 



The Lagrangian density of our (2+l)-dimensional model is given by 



C = + d»od»ä - ^(0 2 - 4f - ^(\a\ 2 - V ^ - p^\a\ 2 + ^Vt (2-1) 

where and a are real and complex scalar fields respectively, with rj^, T] a , \<f,, X a , (3 all real 
positive constants. This Lagrangian density can be obtained from the global version of 
Witten's Z7(l) x U(l) model [15] by restricting one of the complex scalar fields to be real. 

The theory has a global Z 2 x U(l) symmetry and the parameters of the model can be 
arranged so that in the vacuum the Z 2 symmetry is broken, = ±77^ 7^ 0, while the U(l) 
symmetry remains unbroken, |er| = 0. For this symmetry breaking pattern there exist kink 
strings constructed from the <fi field. If the infinite kink string lies along the ?/-axis, then it 
is given by the solution 



= ^tanh(^ V 2 j, a = 0. (2.2) 

The situation of interest is when a condensate of the a field carrying current and chargé 
forms in the core of the kink string. For the infinite string given above, such a condensate 
field takes the form 

a = é {u)t+ky) \a\, (2.3) 

where |cr| is a function of x only with |er| — > as |x| ^ 00. The constant k describes the rate 
of twisting of the condensate along the string, though in the literatúre this is referred to as 
winding, rather than twisting, so we will stick to this common convention. 

A non-zero value of u induces a chargé Q associated with the global U(l) symmetry, and 
the winding k generates a current along the string. It is easy to see that chargé and current 
will have opposite effects on the string, so it is useful to introduce the combination 

X = uj 2 - k 2 . (2.4) 

In the literatúre solutions with x — are termed chiral, whereas solutions with x > are 
referred to as electric and those with x < are called magnetic [H [TU] • 

It is also convenient to introduce the quadratic coefficients for both fields by the definitions 

2 _ Afli 2 2 _ A(7 2 i ŕn r-\ 

m 4, = y v ^ m<7 = Y VcT x ' ^ ' 

Taking into account the form of the condensate (12.31) . a simple analysis of the potential term 
in (12.11) . reveals that the condition for the Z 2 symmetry to be broken in the vacuum is 

f > f- < 2 - 6 > 

Similarly, for the U(l) symmetry to remain unbroken in the vacuum requires that the effective 
mass term (including a contribution from the interaction) for the a field is positive, which 
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results in the constraint 

P > ^4- (2.7) 

Note that both conditions (I2.6P and (12.71) are identical to those that arise in the cosmic 
vorton case [TU] . 

For a condensate to form in the string core we require that the solution (12. 2p with no 
condensate is an unstable solution of the field equations 

<W + ^ - ^ + P<P\a\ 2 = 0, (2.8) 

d^v + ^{\a\* - rfta + Pťa = 0. (2.9) 

Linearizing these equations around the solution (12.21) with a o field of the form 

a = e int e iiujt+ky) \a\, (2.10) 
where again \a\ is a function of x only, yields the eigenvalue equation 



- \a\" + (/3^tanh 2 (^-J - - x)W\ = íŕ\a\ (2.11) 

This is a classic example of a stationary Schrôdinger problém where the spectrum is known 
explicitly [9]. Requiring a negative mode, Q 2 < 0, gives the condition 

A m^(2m^ + m^) 
/?< 2m^ • 

This is similar to the condition derived for cosmic vortons [71 [9], but in that case some 
approximations must be made, whereas here this is an exact result. 

The three conditions (I2.6p .( l2ľ7l) . (l2.12p provide constraints on the constants of the theory 
that must be satisfied for some range of the parameter \. 

3 Exact solutions 

For an infinite státie string, where the condensate field has the form (12.31) . the field equations 
(12. 8p and (12. 9p reduce to the following nonlinear coupled ordinary differential equations 

<P" = <P^f(<P 2 ~ Vl) + PW\ 2 ) (3.1) 

W\" = W\(-x + Y^\ 2 -rt) + P<P 2 )- (3-2) 

If x — then an exact solution to the above equations can be found if [T2| [8] 

H)H! + f- 2 k (3 - 3) 
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For the purposes of the current investigation it is useful if exact solutions are available for a 
range of values of the parameter \. This can only be achieved for the choice 



2/3 = A^ = X a = X. (3.4) 

With this relation between the constants of the theory a simplification occurs and the two 
conditions (I2.6P and (12.71) become identical and reduce to the condition > rrí^. Fur- 
thermore, the final condition (12.121) simplifies to m| < 2m 2 a) so the whole set of conditions 
reduces simply to the constraint 

1 . . 

-^ m 4> < r m a < m f (3.5) 

The exact solutions to (13.1 j) and (13.21) are given by 



(j) = T]^ tanh (x^Jml - m£j , \a\ = \J — ^ ^-sech (x^Jm^ - n\?\ (3.6) 

which are valid for the entire range of x satisfying the condition (13.51) . 

At this stage there are three remaining constants of the theory, that is, rj a , X. Most 
of the analysis which follows can be extended to arbitrary values of these constants (provid- 
ing they satisfy (13. 5p for some range of x) but in order to simplify the presentation and to 
compare with later numerical simulations we shall now fix some specific values. In the penul- 
timate section we shall discuss how other choices for the constants of the theory inffuence 
the results. 

In the chiral case, x — 0, the condition (13.51) reduces to 

To aim for the most generic behaviour possible we set r\ a = \^3r] ( j ) /2, so that the value chosen 
for the ratio which appears in (13.71) is at the midpoint of the allowed range. Finally, the 
overall scales are normalized by setting t]^ = 1 and A = 2. In summary, the constants selected 
in the remainder of this paper (results for more generál choices are presented in appendix 
A) are 

tto = l, ^ = ^3/2, A^ = 2, A ff = 2, (3 = 1. (3.8) 
With the above choice the constraint ( 13.5P reduces to 

-\<X<\ (3.9) 

hence a reasonable range of electric and magnetic solutions are possible in addition to the 
chiral case. The infinite string solutions (I3.6P simplify to 



tanh 1*^7^1 H = , iiíísech í Í^ZM). (3.10) 
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From these exact solutions the result of adding current (x < 0) an d chargé (x > 0) to the 
condensate can easily be seen. Current clearly quenches the condensate, with the maximum 
amplitúde decreasing as x decreases, until finally the condensate disappears in the limit as 
X — > —1/4, which is the lower edge of the allowed range ŕ)3.9p . There is also a moderate 
decrease of the width of the condensate (and kink) as the current is increased. On the 
other hand, chargé enhances the condensate, with the maximum amplitúde increasing as x 
increases. Moreover, the width of the condensate (and kink) increases dramatically as x 
increases, and becomes delocalized in the limit x ~ * 1/4, which is the upper edge of the 
allowed range (13.91) . Similar quenching and anti-quenching effects are seen in the numerical 
solutions of infinite cosmic strings carrying a condensate in the (3+l)-dimensional model 

M- 



4 Kinky vortons 

We now turn our attention to kinky vortons, and apply some approximations, in conjunction 
with the exact results of the previous section, to study circular loops of kink strings carrying 
a condensate. In the subsequent discussion we will use the terms vortons and cosmic vortons 
to refer to kinky vortons and their three-dimensional analog respectively. 
In polar coordinates the ansatz for a circular loop reads 

</> = <f>(r), o = é^ +m) \a\(r), (4.1) 

with the boundary conditions 0'(O) = 0, 0(oo) = 1 and |c|(0) = 0, |cr | (oo) = 0. This 
describes a vorton of rádius R where 4>(R) = 0. The integer N is referred to as the winding 
number. Note that to compare with the straight string the winding rate k depends upon 
both iV and the loop rádius R since k = N/ R. 

The conserved chargé Q of such a field is given by 



1 



Q = — I (áä - aa) d 2 x = 2ttuj I \a\ 2 rdr. (4.2) 
2i J Jo 

The main aim is to determine the properties of a vorton, and in particular its rádius R, as 
a function of the conserved quantities Q and N. 

If the rádius R is much larger than the width of the condensate (and kink) then a section 
through the circular loop can be well-approximated by a section through the infinite straight 
string. This means that approximations of the following form can be used to evaluate a 
number of required integrals 



\a\ p rdr ta KĽ p = R 



l + 4x , fxy/1 - 4x 
seen 



v 



dx, (4.3) 



where the infinite string solution ŕl 3 . 1 j) has been used to determine the cross-sectional con- 
tribution. In particular, we shall later require the exact results 

^ 2(1 + 4 X ) ^ 2(1 + 4 X ) 2 



VT^F 3v^4^ 



(4.4) 
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In the analysis of cosmic vortons analogous quantities to S 2 and S 4 are also required, but 
in that case they can only be computed numerically; although analytical estimates have 
been given [TU] which have a similar form to the exact results Íj4.4p . The main qualitative 
difference is the appearance of the square root in the denominators in (I4.4p . which can be 
attributed to the reduction in the number of spatial dimensions. 
Substituting the above result for Ľ 2 into (14. 2 p yields 

Q = ^u;R^±^=. (4.5) 

In the chiral case, x — 0> then u — k — N/ R, thus (14. 5p reveals that for chiral vortons the 
relationship between the conserved quantities is Q = 4ttN, and is independent of R. Thus 
loops which are initially chiral remain chiral even if the rádius changes. 

The first goal is to determine the rádius R at which the vorton energy is minimized for 
fixed values of the conserved quantities Q and N. One difficulty is that the cross-sectional 
properties of the string (for example, the amplitúde and width of the condensate) vary with 
the vorton rádius. As seen above, this variation is most easily quantified as a function of 
X, and is highly implicit as a function of Q and N. This is a serious obstacle for cosmic 
vortons, but the fact that we have exact expressions available makes the problém much more 
tractable. 

Taking the square of equation (14.5 p and making the substitution uj 2 = x + N 2 /R 2 gener- 
ates the cubic 

16X 3 + 8(l + x + (íŕ _ ^ 1 = o. (4.6) 



R 2 J \ R 2 47r 2 ÄV V 1Q7T 2 JR 2 

Given the physical properties of the vorton, that is the conserved quantities Q and N together 
with the rádius R, this cubic determines x ( as the root which lies in the interval (—1/4, 1/4)). 

The cubic (14.61) allows the variation of x with the rádius R to be easily studied. First of 
all, the cubic confirms that the chiral case is a fixed point, that is, if Q = 4irN then x — 
for all R. Also, as R — > oo the final term in (14 .6p tends to zero, thus x ~ > Oj demonstrating 
that if a vorton is created with a sufiTciently large rádius it will be arbitrarily close to chiral. 
Early studies of cosmic vortons suggested that the chiral limit might be an attractor [BJ but 
more recent investigations [10J have demonstrated that the chiral limit is a repulsive fixed 
point not an attractor. For kinky vortons this can easily be seen from (14.61) as follows. The 
electric regime is given by Q > AirN, in which case the final term in (14.61) is negative and 
X > 0. For large enough R then x ~ 0, but as R decreases the modulus of the final term 
in (j4.fi p increases and therefore so does Xi moving away from the chiral limit. Similarly, if 
Q < AttN, then all coefficients in the cubic (14. 6j) are positive thus x < an d this is the 
magnetic regime. As R decreases then again the modulus of the final term in (14 .6p increases 
and therefore x becomes more negative, again moving away from the chiral limit. This is 
illustrated graphically in Figúre [TJ where curves are plotted showing the value of x as R 
decreases for Q = 1500 and N = 250, 200, 150, 100, 50, 10 (left to right). In each case as R 
decreases the curve moves away from the chiral line x = 0- Each of these curves terminates 
at the energy minimizing rádius, which we now discuss. 
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x 

Figúre 1: The variation of x as t ne vorton loop rádius R decreases for Q = 1500 and 
N = 250,200,150,100,50,10 (left to right). Note that as R decreases each curve moves 
away from the chiral line x = 0- Each curve terminates at the energy minimizing rádius. 



Using the cross-sectional formulae (14.41) and similar expressions for the cross-sectional 
integrals of powers of 0, cf>' and \a\' the vorton energy is found to be 



E 



R 



5% + 4x z 



+ 



R 



iV 2 (l + 4 X ) + 



Atx) 



(l-4x) 



;i+4x) 



(4.7) 



For given values of the conserved quantities Q and N, this defines the energy as a function 
of the vorton rádius R, where x is again given by (14. 6 p . Using this formula the value of R at 
which this energy is minimized can be calculated. This is plotted in Figúre [2] for 1 < N < 300 
and the three values Q = 2000 (dashed curve), Q = 1500 (solid curve) and Q = 500 (dotted 
curve). 

The energy formula (14. 7p together with the cubic (14.61) can be used to show that if either 
Q = or iV = then the energy minimizing rádius is R — 0. Thus vortons require both 
current and chargé, in agreement with the results for cosmic vortons [TU] . 

The above analysis for the vorton rádius as a function of Q and ÍV can be compared with 
the results of numerical simulations of the full nonlinear field theory. Restricting to radially 
symmetric fields of the form (14.1 p the energy can be written as 



E = 2tt 



i:\ 


/2 + M 







D 2 + M 2 



a N 1, .2 3 
r 2 2 4 



- ) > r dr + 



Q 2 



oo 




a\ 2 r dr 

(4.s; 
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Figúre 2: The energy minimizing rádius R as a function of N for the three values Q = 2000 
(dashed curve), Q = 1500 (solid curve), and Q = 500 (dotted curve). The results of a field 
theory energy minimization are presented as filled circles for Q = 1500 and open circles for 
Q = 2000. We note that the solution with Q = 1500, iV = 10 and R = 50 is the smallest 
rádius solution we have been able to find. Vorton solutions will only exist for points on each 
of these curves down to a minimal value of R due to the splitting instability. 
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For fixed values of N and Q this energy has been minimized using a simulated annealing 
algorithm. The resulting loop rádius for the four values N = 10, 100, 200, 300 are presented 
in Figúre [2] as the filled circles for Q = 1500 and the open circles for Q = 2000. From 
these, and other similar results, it can be seen that the analytic calculation compares well 
with the numerical results, when a vorton solution exists. However, if Q and N are not 
sufficiently large then it appears that no vorton solution exists. For example, for Q = 500 
(corresponding to the dotted curve in Figúre [2]) no solutions were found for any value of N 
and for Q = 1500 no solution was found with N < 10. This is because of a splitting instability 
where the effective potential trapping the condensate on the kink loop is not strong enough 
to prevent a separation in which the condensate radiates away to infinity leaving behind a 
bare kink loop that collapses. This instability can not be seen in the analytic calculation, 
which assumes that the condensate is trapped on the kink loop, but is clearly evident in the 
numerical energy minimization process. In principle, it should be possible to numerically 
map out the región of the (Q, iV)-plane in which vortons exist. We have not performed 
this computation due to the fact that it requires significant calculation to determine the 
boundary of instability since it can take a long time to develop when close to the critical 
value. 

5 Field theory dynamics 

In this section we describe the results of numerical simulations of the dynamical field theory 
equations ( 12. 8 j) and ( 12.91) . The numerical algorithm is an explicit finite difference schéme 
which is second order accurate in both the space and time derivatives. For most simulations 
the lattice spacing is taken to be Ax = 0.5, although for some simulations this is reduced 
to Ax = 0.25 when greater resolution is required. The time step is Aí = 0.1 for all the 
simulations presented. Neumann boundary conditions are applied and the number of grid 
points ranges from 401 x 401 to 2801 x 2801 depending upon the size of the vorton that is 
being simulated. 

To study the stability of a vorton to axially symmetric perturbations we wish to consider 
its evolution from an initial state in which the rádius differs from that of the optimal value. 
Of course, the initial field must have the correct chargé Q and winding iV of the vorton we 
wish to study. To create such an initial condition for a circular loop with an initial rádius Rq, 
we use the radially symmetric ansatz (14.11) with the required value of N and where the fields 
through a section of the loop are approximated by the infinite straight string fields (I3.10p . 
Let xo be the value of x used in the initial conditions, which is not that of the originál vorton 
to be investigated, but is obtained by substituting the required values of Q and N, together 
with the initial rádius Rq, into the cubic (14.61) . The initial frequency u Q is then obtained as 

uj = \Jxo + (N/Ro) 2 , where we have made use of the formula (12.41) . 

Vortons with a small loop rádius require less computational resources to study than larger 
vortons, simply because fewer grid points are required to fit the vorton inside the simulation 
región. From Figúre |U and recalling the earlier comments that large values of Q appear to 
be necessary to avoid the splitting instability, it can be seen that the strongly electric regime 
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Figúre 3: The rádius of an electric vorton for < t < 50000. The chargé and winding 
are Q = 1500 and iV = 10, corresponding to an optimal rádius of R = 50, but the initial 
condition has R = R = 70. The inset is a blow-up of this plot for early times < t < 2000. 
It can be seen that the rádius oscillates around the optimal value but the amplitúde decreases 
extremely slowly. 




ABC 



Figúre 4: An electric vorton at t = 0, 8000, 20000. Dark (red) regions are where <fi 2 < 0.2. 
The chargé and winding number are Q = 1500 and = 10, corresponding to an optimal 
rádius of R = 50. The initial condition is a circular loop of rádius R = 50 which has been 
squashed by 10% along the y-axis and stretched by the same factor along the x-axis. 
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44 O 5000 10000 15000 20000 

t 

Figúre 5: The rádius along the positive x-axis of an electric vorton for < t < 20000. 
The chargé and winding are Q = 1500 and N = 10, corresponding to an optimal rádius of 
R = 50, but the initial condition is an ellipse. The inset is a blow-up of this plot for early 
times < t < 2000. Long-lived ringing modes are observed once again. 



is the easiest to study. In particular, there is a solution with Q = 1500 and N = 10, with 
a corresponding value of % = 0.193 and a rádius R = 50, which is the smallest vorton that 
appears to exist in this theory. Field theory simulations of this vorton using the correct initial 
rádius Rq = 50 confirm that the vorton rádius remains constant throughout the simulation. 
To test stability we consider an initial condition in which the rádius is R$ = 70 and therefore 
is substantially larger than the correct value for the given chargé and winding. The evolution 
of the rádius with time is displayed in Figúre [31 The inset shows the evolution of the rádius 
for earlier times < t < 2000, when only a few oscillations of the rádius have taken plače. 
This confirms that the rádius oscillates around the optimal value R = 50 but reveals that 
the amplitúde of the oscillation decreases extremely slowly due to radiation. The full plot 
extends to very large times < t < 50000 and reveals the decaying amplitúde and approach 
to the optimal rádius, but even after over a hundred oscillations the amplitúde is still a 
reasonable fraction of the initial perturbation. The ringing modes are surprisingly long lived 
and the simulations require a large number of time steps (in this example half a million) to 
observe the approach to equilibrium. 

We have performed a number of simulations similar to the one described above, and all 
suggest that when vorton solutions exist they are stable to axially symmetric perturbations. 
In order to test stability to generál perturbations we require an initial condition which 
breaks the axial symmetry, but preserves the chargé Q of the unpertubed solution. Such 
an initial condition is obtained by taking a vorton solution and performing the combined 
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scaling transformation x \— > xfi and y i— > which clearly leaves the chargé Q invariant 
and converts the circular vorton into an ellipse. In the study of cosmic vortons with a 
modified interaction term, elliptic initial conditions were also considered, but in that case it 
appears that these were required to be fine-tuned to ensure an initially homogeneous current 
distribution. As we shall see, our vortons are more robust and do not require any fine-tuning 
of the initial elliptic loop. 

As an illustration, we again consider the electric vorton with parameters Q = 1500 and 
N = 10 yielding the rádius R = 50. The initial condition has squashing parameter fi — 1.1 
to produce the ellipse displayed in Figúre H]A Dark regions are where 2 < 0.2, and hence 
indicate the position of the kink string. The vorton evolves as a rotating and oscillating 
ellipse, with the initial elliptical shape still apparent even at the much later time t = 20000 
Figúre HD. It is expected that the eccentricity of the initial loop decreases with time so 
that the unperturbed circular vortex is eventually recovered. However, as with the axially 
symmetric perturbation, there are extremely long-lived ringing modes. In Figúre [5] we plot 
the rádius along the positive x-axis as a function of time, to demonstrate the slow approach 
towards equilibrium. 

The field theory dynamics described so far has been concerned with strongly electric vor- 
tons, but similar results have also been obtained for chiral vortons. As an example, Figúre O 
displays the evolution (t = 0, 5000, 10000) of a squashed chiral vorton. The unperturbed vor- 
ton has Q = 1500 and N = 119 giving a rádius R = 185 and x ~ 0. The elliptic perturbation 
is given by a squashing parameter fi = 1.1. The chiral vorton is substantially larger than the 
strongly electric vorton studied above with the same chargé Q. Recall that as x decreases 
the kink and condensate width both decrease, so that the ratio of the vorton rádius to kink 
width is significantly larger in the chiral case than in the strongly electric regime. Chiral 
vortons are certainly well inside the thin ring limit. This is illustrated in Figúre [6] where the 
dark región is where (f 2 < 0.5. Note that this value of 0.5 is larger than the value 0.2 used in 
the plots of the electric vorton, and therefore makes the chiral vorton appear thicker than it 
would have done if the same value as earlier had been used. However, for the chiral vorton 
a much bigger grid is required and the results in Figúre M correspond to a grid containing 
1001 x 1001 grid points, but with the same lattice spacing Ax = 0.5. These results confirm 
the stability of chiral vortons and also demonstrate that significant computing resources are 
required even to study kinky vortons in the chiral regime and that similar studies of cosmic 
vortons would be prohibitive. This helps to explain why numerical studies of cosmic vortons 
are currently limited and the only solutions found to date are in a theory with a modified 
interaction term designed to reduce these difficulties. We shall make some further remarks 
on this topic later. 

Analytic and numerical work [TU] on infinite cosmic strings carrying current and chargé 
has shown that if the string is strongly magnetic (x <C 0) then it will develop a pinching 
instability, in which the condensate field is driven to zero at some point on the string and the 
field partially unwinds to produce a string with a lower winding number N. Based on this 
result it has been conjectured that a magnetic vorton with \x\ > Xc will excite the unstable 
pinching mode, reduce its winding number and recover as a vorton which is less magnetic, 
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ABC 

Figúre 6: A chiral vorton at t = 0, 5000, 10000. Dark (red) regions are where (f) 2 < 0.5. 
The chargé and winding number are Q = 1500 and N = 119, corresponding to an optimal 
rádius of R = 185. The initial condition is a circular loop of rádius R = 185 which has been 
squashed by 10% along the y-axis and stretched by the same factor along the x-axis. 



or possibly even electric. We shall now investigate the dynamics of magnetic kinky vortons 
and demonstrate this process has a significant impact on the dynamics. 

An analysis of inŕmite kink strings to a pinching instability is presented in appendix B 
and suggests that such an instability will be present if x < —0.1. This should provide a 
reasonable estimate for the value at which such an instability emerges for kinky vortons, 
but is not expected to be exact as it neglects the curvature of the vorton loop. In fact the 
numerical simulations described below reveal an instability when x ~ —0.08, which is in 
reasonable agreement with the analytic calculation. 

A magnetic vorton will have a larger rádius and a smaller width than its chiral counterpart 
so even larger grid sizes must be used to simulate the magnetic regime. The following 
simulations of magnetic vortons used a grid containing 2801 x 2801 points and a reduced 
lattice spacing of Ax = 0.25, with 80000 time steps per simulation. Once again it is worth 
pointing out that analogous simulations in (3+l)-dimensions would require very considerable 
computing resources. 

First of all we shall consider a magnetic vorton which is expected to be stable. If Q = 1500 
and N = 159 then R = 213 and x = ~ 0.04, which should be within the stable regime. The 
dashed curve in Figúre [TJ shows that the rádius oscillates with a small amplitúde around the 
initial value, suggesting that this magnetic vorton is indeed stable. 

For the values Q = 1500 and N = 214 the vorton rádius is R — 245 and x = —0.08, 
which is expected to be close to the onset of the pinching instability. The solid curve in 
Figúre [7] displays the rádius as a function of time for this vorton. It can be seen that the 
rádius decreases significantly towards the štart of the evolution and subsequently oscillates 
around a much lower value. This is precisely the behaviour expected of an unstable magnetic 
vorton which partially unwinds and recovers to a vorton with a lower value of N. This 
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Figúre 7: The rádius as a function of time for two initially magnetic vortons. The dashed 
curve corresponds to a vorton with initial parameters Q = 1500 and N = 159 giving R = 213 
and x — —0.04. This vorton is stable and the rádius oscillates around the initial value with 
a small amplitúde. The solid curve is for a vorton with initial parameters Q = 1500 and 
N = 214 giving R = 245 and x — —0.08. This vorton is unstable and the rádius decreases 
significantly near the štart of the evolution and then oscillates around the lower value of 
R = 167. 



15 



A 



P, 



C 



Figúre 8: The dark (red) región is where 5ŕ(cr) > 0.05, plotted at the times t = 0, 1000, 8000. 
Initially it can be seen that the winding number is N = 214 but this has reduced to N = 98 
in the two subsequent plots. 



interpretation is confirmed in Figúre [8] where we plot the dark región where 5ŕ(cr) > 0.05 
at times t = 0, 1000, 8000. The advantage of displaying the quantity 5ŕ(cr) is that not only 
can the location of the condensate (and hence the vorton rádius) be observed but so can the 
winding number N, since it corresponds to the number of dashes contained in the vorton loop. 
Counting the number of dashes in Figúre [8]A. confirms that initially iV = 214. Counting the 
number of dashes in Figúre [HJ3 and Figúre [8]C reveals that each configuration has a winding 
number N = 98. A vorton with Q = 1500 and iV = 98 has a rádius R = 167 and this is 
consistent with the rádius around which the solid curve in Figúre [7] oscillates at later times. 
This vorton has x = 0.03 and is therefore electric. These results confirm that the initial 
magnetic vorton is unstable and partially unwinds to convert to a stable electric vorton. 
Note that the pinching instability can only be excited if the axial symmetry of the vorton is 
broken but this is provided by the small perturbation produced by the numerical grid and 
in particular the boundary. This can be seen by the square deformation of the vorton in 
Figúre Ep. 

6 Parameter values and electric instability 

The field theory simulations discussed in the previous section illustrate that the computa- 
tional resources required to study strongly electric vortons are not as substantial as for chiral 
or magnetic vortons. This is relevant for the study of cosmic vortons, where simulations in 
(3+l)-dimensions are close to the limit of current feasibility. It is more likely that cosmic 
vortons can be constructed numerically if they exist in the strongly electric regime. In this 
section we shall discuss how the existence of strongly electric kinky vortons can depend 
crucially on the specific parameter values of the model. 
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Recall from earlier, that vortons do not exist if they are too electric because of a splitting 
instability, where the effective potential trapping the condensate on the kink loop is not 
strong enough to prevent a separation in which the condensate radiates away to infinity. 
The only current numerical construction [TU] of cosmic vortons is in a modified model where 
the interaction term |0| 2 |<r| 2 is replaced by a non-renormalizable interaction term |0| 6 |<r| 2 . 
The effect of this modification is to increase the width of the effective trapping potential and 
prevent the separation of the condensate and the string, thereby allowing solutions which 
are more electric (and therefore of a smaller rádius) than in the originál theory. Below we 
demonstrate that the effective trapping potential, and its dependence on x, can vary with 
the choice of the parameters in the Lagrangian. In particular, the parameters used in this 
paper (chosen to allow explicit exact inŕmite string solutions carrying current and chargé) 
produce a trapping potential which does not decrease as x increases, in contrast to some 
other choices. This suggests that the parameter values used have allowed the construction 
of strongly electric vortons that may not exist for some other parameter values. 

To study the trapping potential we introduce the normalized interaction 

f°° ó 2 \cr\ 2 dr 

Y (y) = J -°°^ 11 (6 1) 

which we refer to as the splitting parameter. In this deŕmition <fi and o are the inŕmite 
straight kink and condensate fields respectively, for a given value of %. By construction we 
have that < Y < 1 and the larger the value of Y the weaker the trapping potential, since 
the condensate samples less of the kink field away from its vacuum value. This is why we 
refer to Y as the splitting parameter, since larger values of Y suggest a greater tendency to 
instability via splitting. 

For the parameter values (I3.8P used in this paper, the exact solution can be used to 
compute T (x) explicitly, yielding the result T (x) = 1/3, independent of x- This is a special 
situation and is not generic. As an example, consider the parameter values used in the study 
of cosmic vortons in [10], namely 

13 3 
V<f> = 1, rja = 2> -V =2' A<J = 10 ' ^ = 2" ( 6 ' 2 ) 

For this parameter set there are no exact solutions available and hence the fields and V (x) 
must be computed numerically. The result is displayed as the solid curve in Figúre [9j for 
the allowed range of x (f° r comparison the dashed line is a numerical computation for the 
parameter set (I3.8P and is a test on the numerical accuracy since the exact result T (x) = 1/3 
is known in this case). Figúre M reveals that T (x) increases with x f° r the parameters (16. 2p . 
and suggests that the splitting instability is more severe for this parameter set, hence making 
it more difficult to find strongly electric (and therefore small rádius) vortons. We have been 
able to compute some kinky vortons with this parameter set but the solutions indeed appear 
to exist only in a regime which is much closer to the chiral limit. This fact, together with 
the absence of exact solutions, makes it a more difficult problém to study in detail. 

As mentioned above, one approach to partially overcome the splitting instability is to 
modify the power in the interaction term [TU]. However, the above results suggest that an 
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Figúre 9: The splitting parameter T (x) f° r three choices of the parameters 77^, r/ a , A^, A CT , j3 : 
solid curve for parameters (16.21) ; dashed Hne for parameters (13.81) ; dotted line for parameters 

(ED. 

alternatíve possibility is to find a more favourable región of parameter space. The parameter 
set ( 13.81) appears to be a better choice than (16.21) but there is no reason to believe that the 
parameters which allow exact solutions are the best choice from this point of view. As a 
further example, consider the parameter set 

r),), = 1, ľ] a = 1, A^ = 3, X a = 2, (3 = 2. (6.3) 

The associated splitting parameter V (x) is shown as the dotted curve in Figúre [Hl for the 
allowed range of x- m this case T(x) decreases as x increases, suggesting that the param- 
eters ( 16.31) could be an improvement over the values (I3.8P in searching for small electric 
vortons. This issue, and particularly its implications for cosmic vortons, is currently under 
investigation. 

7 Conclusion 

In this paper we have introduced and studied kinky vortons, which are (2+l)-dimensional 
analogues of cosmic vortons. The lower dimensional systém allows some exact results to be 
obtained and greatly simplifies both the analytic and numerical treatment. We find that 
there are remarkable similarities between kinky and cosmic vortons, and that many of the 
expected properties of cosmic vortons can be realized in the kinky vorton systém. 
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Our study of kinky vortons has revealed some of the difficulties that are responsible for 
the current limited results on cosmic vortons and has suggested some possible approaches 
to tackle these issues. Investigations of cosmic vortons are currently in progress which make 
use of the results presented here. 

Most studies of cosmic vortons apply a string approximation (see |3j and references 
therein) which should be valid in the thin ring limit, where the vorton rádius is much larger 
than the cosmic string width. The simulations of kinky vortons suggests that this will 
provide a good approximation providing the vorton indeed exists and is stable. However, 
the existence and stability appear to require a detailed analysis of the full field theory 
equations, so some caution may need to be applied when using results derived within the 
string approximation. 

It is perhaps worth pointing out that vorton-like objects occur as classical solutions in 
several field theories describing condensed matter systems [TJ EJ Q2] , so in the future it may 
be possible to study these objects experimentally within this setting. 
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Note Added 

During the final stages of the preparation of this paper an interesting preprint [TT] appeared 
on the arXiv containing some numerical results on the construction of cosmic vortons. How- 
ever, the cosmic vortons in [11] are in a very different regime to those discussed in this paper. 
In particular, the vorton rádius is almost equal to the cosmic string width, so they are far 
from the thin ring limit. Indeed, only very low winding numbers, N < 5, could be computed 
and it was acknowledged that the construction of vortons in the thin ring limit remains a 
numerical challenge. The vortons in [H] were constructed by perturbing away from the 
sigma model limit A^ = A CT — > oo and (3 — > oo and it was noted that such vorton solutions 
are essentially the Skyrmion solutions obtained by the current authors [1] in a related model 
of a two-component Bose-Einstein condensate. 

Appendix A 

In this appendix we provide some formulae for kinky vortons for generál values of the pa- 
rameters A, 7)$ and r/ a , while maintaining the condition A = = A CT = 2(3. If we define 
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(Ver / V<t>) 2 an d X = x/ m then the constraint (J3ľ5]) becomes 



-- a <x<\ - a. 



(A.l) 



and the solutions (13.61) can be written as 



rjtf, tanh 



rriéxJl - a- x 



a\ = VŽrj^i a — - + xsech m^x^Jl — a — x 



(A.2) 



Using these formulae one can show that 

IQttRu (a — | + x) 



Q 



and x satisfies the cubic 

1\ 16N 2 



16 X Ó + 



32 a 



R 2 



X 



A V 1 - ä - x 



, ( M,-i) 2 + 32f«-i)^ + f^y 

2y V 2/^2 ^ 47rÄ y 



(A.3) 



+ 



16(«-Í)V + (a-l) 



— = 
i? 2 



(A.4) 



where i? = m^R. The energy is then given by 



E 



2^ 2 



yl-a-x V 3 



R 



2(2a+l)(2-2a- X )+4 x s 



1 

R 



2 n 



(A.5) 

For the values r?^, = 1 , A = 2 and a = 3/4, as used in the rest of this paper, these expressions 
reproduce those presented earlier. 



Appendix B 

In this appendix we provide an analysis to estimate the value of x below which a pinching 
instability is expected to exist. This is a lower dimensional analogue of a similar calculation 
performed for cosmic vortons in |10j . 

The instability is studied for a straight inŕmite kink string in the magnetic regime. Using 
the Lorentz invariance of the theory we can restrict to the case u — and therefore x — ~ k 2 . 
The task is to consider the stability of the kink string fields (I3.10p 



tanh 



2V 1 - 4 X 



<TQ 



'l+4x , (Xy/1 - 4 X 
— - — sech 



(b.i; 
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An exact analysis of the stability of this solution would require a numerical approach, so to 
make analytic progress we consider a perturbation of the condensate field of the form 

a = a + \a \v, (B. 2) 

where v (y, t) is a complex function of y and t but is independent of x. There are no exact 
negative modes with the above specified x dependence, as can be seen from the linearized 
equation 

\<T \{j%v - d 2 y v) - u\a \" + H<To|(2|a | 2 - ~ + 0o) + = 0. (B.3) 

However, we can consider integrating this equation over the string cross-section by multiply- 
ing by |o"o| and integrating over x. This produces the equation 

T, 2 {d 2 v - d 2 y v) + z/( X S 2 + £ 4 ) + vX á e 2iky = 0, (B.4) 

and we can now search for negative modes of this equation of the form 

u(y,t) = e At e ik y U (y), (B.5) 

where u (y) is a complex function of y and we require that A 2 > so that A is real with 
a positive value corresponding to the required negative mode. The eigenvalue equation for 
u(y), obtained by substituting ( 1B.5I) into (1B.4I) . is 

- u" - 2%kv! + (u + w)E + A 2 u = 0, (B.6) 

where we have defined S = S4/S2. Note that an equivalent method to derive this equation 
is to substitute the perturbation ( 1B.2I) into the expression for the energy and integrate over 
the string cross-section, then require negative modes for the resulting effective energy. 
Splitting u into real and imaginary parts as u = u\ + iu 2 gives the coupled equations 

- u'{ + 2ku' 2 + (2E + A 2 )^! = 0, -u 2 ' - 2ku[ + A 2 u 2 = 0. (B.7) 

Expanding in terms of fourier modes u\ = a\ cos(py) and u 2 = a 2 sm(jry) yields the equation 

V + 2E + A 2 2kp \/ai 
2kp p 2 + A 2 J \a 2 

A non-zero solution exists only if the determinant vanishes, which gives the condition 

A 2 = -(p 2 + E) ± yj S 2 - A X p 2 . (B.9) 
Requiring a solution with A 2 > produces the constraint 

p 2 + 2S + 4x<0, (B.10) 

and hence 

Using the exact expressions (I4.4p this constraint is equivalent to x < — jô- 

The above analysis has neglected the curvature of the vorton and also integrated over 
the string cross-section, rather than solving the full linearized equations, but it is expected 
to give a reasonable estimate of the critical value of x below which the magnetic instability 
appears. 



0. (B.8) 
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